Scalar radiation from Chameleon-shielded regions 
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I study the profile of the chameleon field around a radially pulsating mass. Focusing on the case 
in which the background (static) chameleon profile exhibits a thin shell, I add small perturbations 
to the source in the form of time-dependent radial pulsations. It is found that the chameleon field 
inherits a time dependence and there is a resultant scalar radiation from the region of the source. 
This has several interesting and potentially testable consequences. 
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Chameleons are scalar fields characterized by a pro- 
file that depends on the local matter density, as a conse- 
quence of their coupling to matter fields. Such a coupling 
naturally arises in string and supergravity theories and is 
often encountered in models addressing the phenomenon 
of cosmic acceleration. For instance, several models of 
dark energy introduce a coupling in the dark sector [1], 
while modified theories of gravity typically introduce a 
nonminimally coupled scalar degree of freedom [2]. In 
all these cases the scalar has a nonlinear potential and 
its effective mass depends on the local matter density, 
i.e. it displays a chameleon behavior [3]. It is precisely 
this nonlinear mechanism that allows the field to have 
a nontrivial dynamics on cosmological scales, while es- 
caping stringent constraints from local tests of gravity 
even if it couples to matter with gravitational strength. 
This screening mechanism has been suggested in [3] and 
further explored in [4-6]. 

As analyzed in [7-9], Birkhoff's theorem ceases to hold 
in these models and there are several interesting conse- 
quences of this fact. An intriguing and yet unexplored 
feature of these theories is that the coupling between 
the chameleon and matter fields allows a radially pul- 
sating mass to directly propagate a disturbance into the 
space surrounding it. Interestingly, this would have sev- 
eral consequences which in principle are testable. On the 
astrophysical side, compact sources could emit scalar ra- 
diation and the rate of energy loss from a binary black 
hole would be bigger than in general relativity. Also, the 
dynamics of core collapse might be significantly modi- 
fied in the presence of a chameleon field. Moreover, in 
conformally coupled theories, the particle masses and the 
coupling constants depend on the chameleon field and, if 
the profile of the field is time-dependent, these quantities 
will all inherit a time dependence. 

These are all interesting scenarios which require a de- 
tailed analysis of the chameleon profile around a time- 
dependent source. In this Letter I specialize to the case 
of an isolated radially pulsating mass. I start with the 
time-independent case and determine the profile of the 
chameleon around a static spherically symmetric mass 
following [3, 5]. Successively, specializing to those cases 
in which the chameleon profile displays a thin shell, I 
add a time-dependent perturbation to the source mass 
and study the corresponding time- and scale-dependent 



perturbation to the static profile. 

Let us consider theories described by the action 



^R-l{d^f-V{^) 



9^ 
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where </> is the chameleon field, ipm are matter fields and 



q(i) 



(2) 



is the metric describing the geodesies for the i-th matter 
field. Pi are dimensionless coupling constants, they can 
be assumed to be of order one and in general they may 
be different for each matter species. 

Varying action (1) with respect to the scalar field one 
obtains its equation of motion 



□0 = y^ 



3i(t)/Mf 



(3) 



For nonrelativistic matter the trace of the energy- 
momentum tensor reduces to the energy density. 



g'^^T/^y = p^*^; we shall use a conserved energy density 



defined as p^*^ = e^^^^^^^ p^^\ Concentrating on a single 
matter component (therefore dropping the i indices), we 
can identify the right-hand side of (3) with the following 
effective potential 



V^^{(P) = ¥{(/)) + pe 



(4) 



i.e. the field profile depends on the local matter distri- 
bution. 

Background profile Let us start from the static case, 
solving for the chameleon profile in the presence of an 
isolated static spherically symmetric mass of density pc 
and radius Rc^ immersed in a background of homoge- 
neous density pc- We shall work in the weak gravity 
regime, assuming that the Newtonian potential is small 
everywhere and the backreaction of the energy density 
in the field (f) is negligible. Therefore the metric can be 
approximated to the one of a Minkowski spacetime, i.e. 
9i^u ~ Vi^iy- This static scenario has been studied in [3-5]. 
The energy density coupled to the chameleon is 



Po{r) 



r <Rc 
r>Rc 



(5) 



We consider inverse power-law potentials of the form 



(6) 
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where M has units of mass and n is a positive constant, 
n > 1. The runaway nature of the potential is important 
for the chameleon mechanism [3] and the inverse power- 
law realization (6) arises in supersymmetric models as 
well as from nonperturbative effects in string theory and 
is desirable in quintessence models [1, 10]. It is also the 
form of a potential characteristic of f{R) theories in the 
Einstein frame. In what follows I use the approximation 
I3(j)/Mp <C 1 which holds through the expansion history 
as shown in [3, 4]. 

For a given value of the background density, p^, the 
mass around the corresponding minimum of the effective 
potential, (pi^ has the form 



ml = n{n + 1) 



Ml 



I will indicate with {(pc^^c) and (^q^^g) the value of 
the field and the mass at the minimum of the effective 
potential respectively for pi = pc and Pi = pc- 

The full problem can be solved numerically, however 
with few approximations one can easily obtain an ana- 
lytical solution [3-5]. For the purpose of this Letter it is 
not relevant to review the details of the solution, but sim- 
ply notice that depending on the parameters of the model 
there are different types of profiles and it will suffice to 
separate them in two classes, as follows. Assuming that 
the field will eventually settle into the external minimum, 
(j)Q^ far enough from the source, the field can remain close 
to <pQ also inside the source or, if the perturbing effect 
of the source is strong enough, the field can be driven to 
the internal minimum 0c inside the source. It depends 
on whether the potential energy cost for not lying at the 
minimum inside the source outweighs the gradient energy 
gain for not introducing a profile in the field. In other 
words, whether the source is strong enough to perturb the 
field profile depends on the relative magnitude of 0g — 
and 6/3Mp$c7 where is the Newtonian potential at the 
surface of the source. If (0g — <l^c)/{^PMp^c) 1, the 
field remains close to (j)Q outside the source, close to the 
minimum (j)c inside the source and interpolates between 
the two minima over a thin radial region ARc close to the 
surface of the mass. This solution is known as the thin- 
shell configuration and more details about its form and 
the conditions under which it is achieved can be found 
in [3-5, 11]; a complete classification of the different con- 
figurations can be found in [4, 5]. 

The thin-shell configuration is the one of interest in 
cosmology since the scalar field can have a nontrivial dy- 
namics on cosmological scales, while the fifth force that 
it mediates on local scales is suppressed by the thin-shell 
coefficient 



0G - 0c 



< 1 



(8) 



Rc 6l3Mp^c 
This way, the chameleon field can easily hide from local 
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FIG. 1: Chameleon profile ip (left) and mass squared (right) in 
presence of an isolated static spherically symmetric source of 
radius Re- The parameters of the model are: n — 1^ pG / pc — 
2 • 10"^ and {mcRc)~^ = 0.01420777 as in [11]. 



tests of gravity while, for instance, sourcing the cosmic 
acceleration on large scales [3]. 

I will use the static profile as background for the time- 
dependent analysis, therefore it is desirable to obtain a 
precise solution by numerical methods, treating the prob- 
lem as a boundary value problem, with boundary condi- 
tions 



dr 



(r = 0) = , 0(r ^ oo) = 



(9) 



I use the relaxation algorithm [12] and rescale the field 
and coordinates defining (f = ^/^c and x = r/Rc. The 
equation for the chameleon reads 



dx'^ 



2 dip 



{rricRcY 



1 



Pq{x) 



1 



,n+l 



(10) 
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where I have used 



<C 1 [3]. For appropriate choices 
of the parameters it is possible to recover the thin-shell 
configuration. I restrict to these cases and show a rep- 
resentative solution for the profile and the mass of the 
chameleon in Fig.l. One can notice a sharp transition 
in the field and mass values over a thin shell; other- 
wise the values are, to good approximation, constant 
with {(p^m) ^ (l,mc) in the inner region {x < 1) and 
{(p^m) ^ {^Gi^g) in the outer region {x > 1). 

Time- dependent profile Let us now add small, time- 
dependent perturbations to the source mass. In particu- 
lar, let us specialize to the case of radial pulsations and 
implement them keeping the mass of the source fixed and 
considering oscillations of its surface around the radius 
Rc. Including the time dependence, the chameleon now 
obeys the following equation 



dx'^ 



2d^ ^ [rricRc? 
X dx n + 1 



p{x) 



1 



,n+l 



(11) 



where r = t/i?c is a dimensionless time. 

We shall decompose the field and the density in a 
background static component and a time- and scale- 
dependent perturbation, respectively (p{r^x) = (po{x) + 
5(p{r,x) and p{r^x) = po{x)[l -\-5{r^x)]] (po{x) satisfies 
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FIG. 2: Time- dependent perturbation to the static chameleon 
profile as a function of radial distance from the source center 
on an equatorial section at time t — 13i?c for the model (16) 




Eq. (10). The perturbation will obey the equation 



{rricRcY 



2dSip 
X dx 



1 



^0 



(12) 



Linearizing in 5(p one gets 

d'^difi d^d(fi 2dd(fi {rricRcf 



dx^ X dx 



Sip - 



n - 



1 



8. 
(13) 

This is an inhomogeneous Klein-Gordon equation with 
a scale-dependent mass given by the mass of the back- 
ground field (pQ. We shall solve it for r > with the 
following initial and boundary conditions 

(5(/:)(0,r) = = StJ(/9(0,x) , lim 5lp{t,x)={). (14) 

As stated above, we hold fixed the mass and geometry 
of the spherical source, and vary its volume via small 
oscillations of the surface. In other words, we assume 
that its radius oscillates around Rc with frequency cjq, i.e. 
Rc[l + Ji?c/^c sin (cjo^c'?")] with SRc/Rc <C 1, therefore 
the density perturbation can be written as 



6{r,x) = 3--^ sm{uJoRcr)0{r)0{l 
Rc 



(15) 



where 0{ii) is the Heaviside step function. 

I have solved numerically both the full (12) and linear 
(13) equation and found that, as long as 5(p)/(p) <C 1, the 
linear solution is a very good approximation to the full 
one. Therefore in the remaining we shall focus on the 
linear equation (13). I show some representative features 
of the solution in Figures 2, 3 and 4 for the model 



n = 1, pq/pc = 2 • 10-^ {mcRcY 
(JoRc = 2, SR/Rc = 0.01 



0.01420777, 
(16) 



where I have used the numerical solution of Eq. (10) for 
the background profile cpo. As can be noticed from Fig- 
ures 2 and 3, the oscillations propagate outside the source 
in the form of ripples in the profile of the chameleon, i.e. 



FIG. 3: Time-dependent perturbation Sp to the static 
chameleon profile as a function of time and radial distance 
from the source center for the model (16). 



there is an outgoing scalar radiation. Looking at Fig. 4 it 
can be noticed that the field responds to the pulsations of 
the source differently in the region inside {x < 1) and out- 
side (x > 1) the source. The amplitude of the oscillations 
in the region outside the source is smaller and the ratio 
of the amplitudes in the inner and outer region depends 
on the parameters of the model, in particular on the ra- 
tio of the masses rric/mQ. Performing a Fourier analysis 
of the oscillations we find a superposition of two modes 
of oscillations, related to the frequency of the source ujq 
and to the internal mass rric. In the right panel of Fig. 4 
the two modes can be noticed in the inner region, while 
in the outer the mode uoq dominates over the other. 

With few approximations we can find an analytical so- 
lution to Eq. (13). It is a Klein-Gordon equation with 
a space-dependent mass which can be well approximated 
by a step function, as evident in Fig.l . Therefore we can 
separate the spatial interval in two regions with constant 
mass, (i.e. m ~ nic for x < 1 and m ~ tuq for x > 1), 
and solve for two standard Klein-Gordon equations. We 
then match the solutions at the boundary x = 1. Let us 
perform a Fourier transform with respect to the time co- 
ordinate t and solve the equation in the frequency space. 
Depending on the value of the frequency cj, we have to 
solve a spherical or modified spherical Bessel equation, 
homogeneous in the outer region and inhomogeneous in 
the inner region. After imposing the boundary condi- 
tions and matching the solutions at x = 1 we find the 
following expression for the Fourier transform of Sep 



5(p ■■ 



Sl 
|AcP 



sinh(Aca:) 

a; 



xcosh(Acx) — 
cosh(A,) - 



-Ag(--i) 



X <1 
X > 1 



^ (17) 

where = ^J{miRcY — (w-Rc)^, i = c, G is real or imag- 
inary depending on the value of w. The source term is 



(jJqRc 



n+1 Rc {oJoRcY - [ojRc? ' 



(18) 
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FIG. 4: Left: Perturbation to the profile of the chameleon as 
a function of the radial distance at time t/Rc = 5. Right: 
Perturbations to the profile of the chameleon as a function of 
time at the distances r/Rc = 1 (red thick curve) and r/Rc = 4 
(black thin curve). Both plots correspond to the model (16). 



In deriving (17) I have approximated the mass with a 
step function. While a more precise solution would model 
also the transition from rric to mo over a finite thin shell, 
this transition was properly taken into account in the 
numerical solution and I have found good agreement be- 
tween (17) and the Fourier transform of the numerical 
solution. We notice that there are two resonances at 
uj = ujQ and uj = rric as found in the numerical solution. 
I derived the equivalent of (17) for the case in which the 
chameleon does not form a thin shell and = Vr. 
The frequency spectrum at a given point in space has 
a different shape and the resonances are now at ujq and 
mo; therefore the resonance at rric signals that the field 
is close to 0c inside the source, i.e. the screening mecha- 
nism is at work. 

The oscillations in the chameleon carry away energy 
from the source. We can estimate the rate of energy loss 
by integrating the flux of the chameleon over a sphere of 
large radius r and dividing by the energy of the source. 
This quantity would be useful to assess the detectabil- 
ity of the radiation, however the result depends signifi- 
cantly on the type of source and model parameters; there- 
fore, for a meaningful estimate one needs to model prop- 
erly physical sources, perhaps extending the formalism 
to the relativistic regime characteristic of strong field 
systems. I intend to address this in future work, how- 
ever a preliminary estimate of the energy loss rate (aver- 
aged over one period of the source) for model (16) gives 
^ O(10~^^)s~^. This value is smaller than known 
rates for binary pulsars [13], although model (16) is not 
strictly comparable to binary pulsars. Once the analysis 
is properly extended to strong field systems, the energy 
loss rate might get closer to the observational threshold, 
offering the opportunity to constrain chameleon theories 
with binary pulsars or core-collapsing objects. 

Conclusions I have analyzed the profile of the 
chameleon field around a radially pulsating source. As 
expected, the pulsations perturb the surrounding space 
inducing ripples in the chameleon field that propagate 
as scalar radiation from the source. The physical con- 
sequences that could give rise to observable features are 



several. For instance, energy could be drained from bi- 
nary black holes faster than in general relativity, the 
ripples in the chameleon profile would induce time- 
variations of masses and fundamental coupling constants 
and the hydrodynamics of compact objects, in particular 
the core collapse, could be significantly modified. Finally, 
it is possible that when the time dependence of the source 
cannot be treated as a small perturbation, the changes 
in the profile of the chameleon would destroy the thin 
shell of the background profile. All these features are 
potentially observable and could offer new ways of test- 
ing chameleon type theories, complementary to existing 
tests [14, 15]. 

In this Letter the analysis was restricted to the case of 
a radially pulsating isolated source, modeled as a homo- 
geneous solid sphere with small oscillations of the sur- 
face and the calculations were carried out in the weak 
field limit. Although simplistic, this scenario is a good 
approximation to many models and allows us to identify 
potentially observable features. The results are promis- 
ing and motivate further work in this direction, with a 
more complex modeling of the sources. In particular it 
would be useful to extend the formalism to relativistic 
cases and to nonspherical configurations. 
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